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The capture and subsequent in-spiral of compact stellar 
remnants by central massive black holes, is one of the more 
interesting likely sources of gravitational radiation detectable 
by LISA. The relevant stellar population includes stellar mass 
black holes, and possibly intermediate mass black holes, gen- 
erally on initially eccentric orbits. Predicted detectable rates 
of capture are highly uncertain, but may be high enough that 
source confusion is an issue. Foreground events with rela- 
tively high signal-to-noise ratio may provide important tests 
of general relativity. I review the rate estimates in the litera- 
ture, and the apparent discrepancy between different authors' 
estimates, and discuss some of the relevant uncertainties and 
physical processes. The white dwarf mergers rate are un- 
certain by a factor of few; the neutron star merger rate is 
completely uncertain and likely to be small; the black hole 
merger rate is likely to be dominant for detectable mergers 
and is uncertain by at least two orders of magnitude, largely 
due to unknown physical conditions and processes. The pri- 
mary difference in rate estimates is due to different initial 
conditions and less directly due to different estimates of key 
physical processes, assumed in different model scenarios for 
in-spiral and capture. 
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I. INTRODUCTION 

The characteristic frequency to which LISA [1,2] is sen- 
sitive is comparable to the orbital frequency at the inner- 
most stable orbit of Schwarzschild black holes of mass 
Mbh ^ IO^Mq and associated Schwarzschild radius, 
rs ■ One of the more promising "guaranteed" sources for 
LISA is the gravitational radiation from the final stages 
of coalescence of low mass (1-1000 Mq) compact objects 
with low mass massive black holes [Mbh ^ 10^^^ Mq). 
White dwarfs, with masses of 0.5-1.3Mq; neutron stars, 
with masses of 1.4Mq; and, stellar mass black holes, with 
masses from order solar masses upwards of IOOMq, may 
all coalesce with central supermassive black holes. The 
stellar population on the most tightly bound, low an- 
gular momentum orbits in the cusp around the central 
supermassive black hole is depleted by prompt swallow- 
ing of the stars "initially" on such orbits. The net flux of 
stars into the black hole is then set by the rate at which 
stars (stellar remnants) can enter the "loss cone" [3] on 



orbits such that gravitational radiation will shrink the 
semi-major axis more rapidly than other processes can 
increase the energy of the object on that orbit. Stellar 
dynamical processes, and possibly star formation pro- 
cesses, determine the instantaneous and mean rate at 
which (compact) stars can enter the gravitational radia- 
tion loss cone and evolve dynamically to coalescence with 
the central black hole. 

Previous papers [4-11] considered the likely rate for de- 
tectable signals from degenerate compact objects in the 
cusps of normal galaxies, coalescing with central black 
holes such as the one inferred to be present in the Milky 
Way. Estimates for detectable signal rates are of the 
order one per year and higher, from these sources, but 
the estimates are sensitive to systematic uncertainties 
in the population contributing to the signal, with very 
large (many orders of magnitude) cumulative formal un- 
certainties in the expected signal rate. Since high rates 
of coalescence are self-limiting (through depletion of low 
mass compact objects at high coalescence rates; or the 
growth of the primary through coalescences, to the point 
where the frequency of the innermost stable orbit is too 
low for LISA to be sensitive to further coalescences) , the 
formal uncertainty in the detectable rate is skewed to 
lower rates of coalescence than the "canonical" estimates 
in the literature. The true rate of detectable coalescences 
is sensitive to the assumptions about the mass function 
of the primary black hole, the number density of stars in 
regions where the appropriate mass black holes are found, 
the mass and number of compact remnants available for 
coalescence over a Hubble time in such regions, and the 
large scale structure of the stellar system in which the 
central supermassive black hole is embedded. 



A. Assumptions 

A number of similar assumptions are made by authors 
in estimating coalescence rates: 

• Galactic cusps: Galaxies tend to have "cuspy" , 
non-isothermal centres. Observations and theory 
suggest that real galaxies have rising density pro- 
files at small radii with: 
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and 7 = 3/2 + p, p e [12-16]. Cosmo- 

logical simulations suggest that a moderate cusp 
forms from initial condition cold collapse in hierar- 
chical models [17], with the cusp possibly modified 
by subsequent black hole growth and interaction 
[15,12]. Alternatively, one may assume that the 
cusp grows with the black hole from an initially 
dense, relaxed cluster of stars, with constant den- 
sity profile at small radii. 

Relaxation: An implicit assumption is that cen- 
tral SMBH are near ubiquitous, with Mbh ^ 
Mcai-Sph oc agg^i [13,18,19], where a is the disper- 
sion of the stellar spheroid, defined at some suitable 
radius far from the black hole. Inside the black hole 
effective radius, rh ~ GMbh/o'^, the velocity dis- 
persion is rising like lim^^o oc r~^^'^ , due sim- 
ply to the Keplerian potential of the central black 
hole. Additional assumptions may be made about 
the stellar velocity anisotropy, /3(r), and the flat- 
tening or possible triaxiality of the stellar density 
profile at small and large radii [12,20,21]. 

Note that the census of low mass supcrmassivc 
black holes Mbh ~ 10"^ - lO'^M© of particular in- 
terest for LISA is not available, the existence and 
properties of this populations is extrapolated from 
the local population (Milky Way and M32), from 
the higher mass SMBH observed, and from circum- 
stantial evidence from low luminosity active galac- 
tic nuclei (AGN), consistent with accretion onto 
lower mass SMBH. 

Stellar population: The relaxation time of the 
stellar population is tr oc a^{r)/p{r), constant for 
p — 0, which thus defines the critical cusp. 

Most relevant galaxies are not relaxed in the cen- 
tre, from observations. However, low mass galaxies 
are more likely to have high stellar densities at rh 
and steep cusps, and therefore may have short re- 
laxation times inside r^. The characteristic galaxy 
has luminosity L* . The galaxies most likely to have 
relaxed cusps with central supcrmassivc black holes 
in the right mass range for LISA are sub-L* galax- 
ies, less luminous than the characteristic luminos- 
ity. It is possible that the stellar population was 
initially relaxed as densities were higher and high 
mass stars were present in galactic centers at early 
times. The current relaxation time may be longer 
than the initial relaxation time, and relaxation pro- 
cesses may have imprinted on the initial conditions. 

Mass function: A general assumption that is 

made is that the stellar population is normal, with 
a Salpeter or Scalo initial mass function. This may 
of course be; wrong, and it has been suggested that 
the initial mass function (IMF) in galactic centres 
or nuclear star forming regions is biased to high 
masses (see Kroupa 2001 for review), which would 



substantially increase the estimated LISA rates. A 
further assumption is made about the white dwarf 
(WD) mass function, as a function of the progeni- 
tor zero-age mains sequence (ZAMS) mass (usually 
assumed to be independent of metallicity); the cut- 
off ZAMS mass for neutron star formation (also 
usually assumed to be independent of metallicity); 
and the cut-off ZAMS mass and mass function of 
stellar mass black holes formed (which ought not 
to be independent of metallicity, see Figer, this 
Proceedings). Binarity may in principle be dif- 
ferent in nuclear star clusters than is observed in 
the field, either due to initial conditions, or due to 
dynamical modifications of the binary population. 
Since mass transfer processes in close binaries af- 
fect the formation channels for compact objects, 
this will affect the inferred rate of coalescence. In 
particular, binary interactions may reopen forma- 
tion channels for classes of compact objects which 
arc not available for single stars, reinjecting a pop- 
ulation that might naively be assumed to be de- 
pleted. Assuming a Salpeter IMF and some upper 
mass truncation, the number fractions of different 
compact object populations are then of the order: 
fwD - 0.1- 0.3, fNS ~ 10-Vwcfc, fsH - 10'^- 
fkick is the fraction of neutron stars receiving small 
enough a kick as not to be displaced from the cen- 
tre of galaxies containing low mass SMBH, estimat- 
ing from the globular cluster population of neutron 
stars, we might infer fkick ^ 0.1. 

• Compact Remnants: Compact remnants can 
then coalesce with SMBH through gravitational ra- 
diation or direct capture. Capture into eccentric 
orbits with pericentres ~ few rs which can "grind 

down" (whereby the orbital semi- major axis shrinks 
gradually) through episodic emissions of gravita- 
tional radiation, before they are upscattered by a 
subsequent stellar encounter, provides the sought 
after periodic signal in the final stages of coales- 
cence (see Sigurdsson & Rees 1997 for discussion). 

• The loss cone: The loss cone - 6ic{r) - is, by 
definition, the region where the time scale for co- 
alescence through gravitational radiation emission 
is shorter than time scale to diffuse out of loss cone 
through random walks or large scatterings in angu- 
lar momentum J (or, rarely, energy, E). 

The rate is sensitive to number density profile 
of stellar population in the center of galaxies, 
p oc r-3/2-P, to the mass function and to any 
anisotropies in the stellar distribution function. 

B. Past calculations 

A number of authors have estimated the rate of co- 
alescences of compact objects into central supermassive 



black holes, using a number of different formalisms and 
assumptions about the underlying "initial conditions" : 

• Formalisms: The primary formalisms used follow 
Prank & Rees (1976), Lightman & Shapiro (1977, 
see also [14,15,26]) and Henon (1973): 

— Hils & Bender (1995) used the Lightman & 
Shapiro formalism calculating the coalescence 
rate in M32 like galaxies (nucleated dwarf el- 
lipticals with high density cores, short core re- 
laxation times and an evolved stellar popula- 
tion). The rate calculated by Hils & Bender 
provides both a useful benchmark for other 
people's calculations, and a good estimate for 
the range of event rates LISA may observe. 
The Hils & Bender rate for white dwarf-SMBH 
coalescence in M32 is 1.8 x 10~* events per 
year. 

— Sigurdsson & Rees 1997, and Sigurdsson 1997, 
1998 used the Frank & Rees formalism, with 
the refinement of estimating also the large an- 
gle scattering rate, which may contribute sig- 
nificantly to the total rate in the "pinhole" 
regime, as opposed to the diffusion regime im- 
plicitly assumed by other formalisms. The 
pinhole regime assumes that change in angu- 
lar momentum of individual stars is large per 
unit dynamical time, compared to the width 
of the loss cone in angular momentum space; 
the diffusion regime applies to the opposite 
limiting case of small changes in angular mo- 
mentum per star per unit crossing time. S Sz 
R '97 found a rate for M32 like galaxies consis- 
tent with Hils and Bender, with a somewhat 
higher net coalescence rate if large angle scat- 
tering is included. S & R found that the cos- 
mological rate is dominated by low mass nu- 
cleated spiral galaxies, like the Milky Way, ex- 
tending to lower (bulge) mass galaxies. Such 
galaxies are more numerous than M32 ana- 
logues [23,24], may have higher central den- 
sities, lower SMBH masses (and surprisingly 
higher coalescence rates) and ongoing star for- 
mation that replenishes the compact object 
population. 

Including large angle scattering, S & R found 
an event rate of 3 x 10~^ per year for white 
dwarf-SMBH coalescence in M32. Excluding 
the large angle coalescence their rate agrees 
with Hils and Bender. 

— Miralda-Escude & Gould (2000) used essen- 
tially the Frank & Rees formalism with the 
assumptions of Bahcall & Wolf. They focused 
on low mass black hole coalescence, with the 
assumption of a relaxed central populations, 
finding a steady coalescence rate for Milky 



Way-like galaxies of ^ 10 ^ per year. In 
comparison, S & R found an "initial" rate 
for LMBH-SMBH coalescence of ~ 10-"^ per 
year in cuspy spirals, declining to a sustained 
merger rate of ~ 10~^— 10~* per year, depend- 
ing on the assumptions about the cusp struc- 
ture and total black hole population. With 
high initial LMBH coalescence rates, the ini- 
tial black hole population is rapidly depicted 
(predicting a burst of coalescence at high- 
moderate redshift) and then settles to an 
asymptotic lower rate as replenishment of the 
inner population matches the steady state de- 
pletion rate. 

— Ivanov (2002) used the Lightman & Shapiro 
formalism, and correcting for the different WD 
mass function, gets a rate lower than but con- 
sistent with Hils & Bender for M32 like galax- 
ies. Note that the statement in that paper 
that S & R get a much lower rate than Hils & 
Bender is wrong and is directly contradicted 
by the S & R paper. Some care must be 
made in comparing the rates since the models 
are scaled to different structural parameters, 
which are not independent. Also, differential 
and integrated rates are of course not directly 
comparable. 

— The current state-of-the-art models arc by 
Freitag (2000,2002), using a time cxphcit 
Monte Carlo realisation of the stellar clus- 
ter surrounding the central black hole, doing 
a Fokker-Planck approximation evolution of 
the stellar dynamics, including SMBH coales- 
cence and growth, following the Henon for- 
malism. The assumptions in Frcitag's work 
are typically for a very massive and compact 
isothermal cluster, with substantial structural 
evolution due to black hole growth and stel- 
lar evolution. As a consequence, the rates in 
Freitag's models tend to be strongly time de- 
pendent, but they are consistent with those 
of other authors, scaling Freitag's initial con- 
ditions to the appropriate late-time quasi- 
steady state values typically approximated by 
other work. 

All the calculations assume some 7 (or p) , scale to 
Mbh and one of cr, or the number density of stars 
at the "radius of influence" of the black hole n*(r/j). 
Inspection of the published papers reveals that the 
dominant source of the difference in the rates is the 
choice of "initial conditions" ; principally the choice 
of evolved or unevolved stellar population (IMF, 
BH mass function) , initial density and whether the 
population is initially relaxed or not. All assume 
isotropy and sphericity. 
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• Normal galaxies: Unrelaxed, nucleated sub-L* 

spirals dominate the rate. 

• Sustainability:Very high coalescence rates are 
self-limiting - they deplete the compact object pop- 
ulation or grow the SMBH to larger mass, moving 
out of the LISA sensitivity range. The coalescence 
rate declines with higher mass SMBH, primarily 
due to the anti-correlation with stellar density and 
correlation with stellar dispersion. 

• Aside - runaway growth: we can check for 
runaway growth under these assumptions, to sec 
whether the growth of the loss cone with SMBH 
mass can trigger rapid mass growth, since adding 
some mass SM ^ Mbh causes the loss cone to ex- 
pand. Runaway growth occurs if the time averaged 
occupancy number of additional stars in the addi- 
tional slice of loss cone, 5d{SM), is greater than 
one. Trivially we find SOic/dM oc 9ic, implying 
exponential growth, but on inspection, with very 
low growth constant, so effectively there is linear 
growth of the SMBH due to the expansion of the 
loss cone, with a low growth rate. 

The exception is for high m^,/MBH where S9{SM) 
can be large and therefore approach runaway con- 
ditions. So, in this approximation, we expect run- 
away growth for the lowest mass seed black hole 
only, which saturates when Mbh 3> , or, if there 
is a substantial population of very massive compact 
objects in the central regions (m* ^ lOOM©). 

C. Predictions 

WD merger rates are uncertain by about a factor of 
three. The primary sources of uncertainty are the rel- 
ative number density of WDs given the observed stellar 
luminosity density; and, whether the inner cusps are typ- 
ically relaxed or not. 

NS merger rates are bounded above to be no more 
than approximately few % of the WD rate, unless nuclear 
IMFs are very skewed to high masses, since the number 
fraction of NS from a normal evolved stellar population is 
small, and the NS are not massive enough to mass segre- 
gate efficiently in the cusp. Since they are more massive 
than typical WDs, they may still contribute significantly 
to the observed rate. However, if NS natal kicks are 
ubiquitous, then the NS population is displaced from the 
inner cusp and the number density in the inner cusp is 
small and the rate negligible at all times. 

BH merger rates are uncertain by at least two orders 
of magnitude and in reality arc likely highly time de- 
pendent. Critical issues is whether the cusp is relaxed, 
allowing an infusion of low mass BH from the outer cusp, 
what the number fraction of low mass BH formed is, and 
whether it varies strongly with the metallicity of the star 
forming gas, and whether star formation persists in the 



inner nuclear region at late times (as seems likely given 
observations of the Milky Way). 

• LISA detection: All predictions in the literature 
are consistent with LISA detections. That is, de- 
spite the range of coalescence rates, even the most 
conservative rate estimates predict LISA will see 
multiple events during its planned operating life- 
time. 

• Source confusion: The highest rate predictions 

imply possible source confusion. 

• The uncertainties in the rate estimates are primar- 
ily assumption driven, dominated by the initial con- 
ditions assumed, not the calculation method. In 
principle, the true initial conditions are constrain- 
able by observations. 

• Other processes: Some significant physics are 

still not included in the models. There is some rea- 
son to believe that the missing physical processes 
will not substantially affect the rate, and will if any- 
thing tend to lead to enhanced coalescence rates, 
with one possible exception. Significant tangential 
anisotropy of the inner stellar cusp could suppress 
the coalescence rate. 

Other issues of concern include: 

Direct capture vs spiral-in. 

Gravitational radiation detection requires orbits 
that spiral in, gradually shrinking the semi-major 
axis over O(IO^) orbits. The capture is into or- 
bits with peribothrons of few rs (e ~ 0.999 — 
0.9999999), but not directly into event horizon, or 
rapidly perturbed into event horizon. Some frac- 
tion of those orbits are scattered down, across the 
event horizon, by stellar encounters near apoboth- 
ron. 

Estimates of direct loss fraction varies substantially 

1/3 or more). There is some concern that this 
fraction could be very large, eliminating a large 
fraction of the estimated sources. Calculations sug- 
gest that the loss fraction is in fact modest [28,5,9]. 

Tidal disruption estimates, similar coalescence pro- 
cesses apply to main-sequence stars, but those stars 
generally undergo tidal disruption far outside the 
Schwarzschild radius (except for the most massive 
SMBH) and we should observe stellar tidal disrup- 
tion events [29,30]. The tidal disruption rate scales 
approximately with the compact object coalescence 
rate (but see Miralda-Escude & Gould 2000), typi- 
cally rcoU for stars is much larger than for compact 
stellar remnants, so the density distribution of stars 
"breaks" at larger radii in the cusp. 

Observationally, the inferred tidal disruption rate 
is consistent with maximum predicted rate Don- 
ley et al (2002). Naively this would imply that 
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the loss cone is generally full. Donley et al also 
find that the inferred tidal disruption rate is signif- 
icantly higher for AGNs, suggesting that dynamical 
processes refilling the loss cone, or nuclear star for- 
mation, may substantially enhance the rate, con- 
tributing strongly to the integrated cosmological 
rate. 



II. FUTURE ISSUES 

A number of open issues remain to be explored in more 
detail. 



A. Brownian wandering 

The formalisms used generally assume the SMBH is 
fixed centrally. 

In reality there can be substantial "Brownian motion" 
of the central SMBH, particularly for the lower mass 
SMBH of particular interest to LISA. In general we have 
some black hole velocity, vbh; mean black hole displace- 
ment from the center of mass of the stellar population, 
{fBH)] and some time scale on which the Brownian mo- 
tion occurs, Ibh- 

If the stellar population is isothermal and relaxed, then 
the problem is trivial, the SMBH is in equipartition and 
moving as a thermally excited massive particle in a har- 
monic oscillator potential with, Viso-SH = ym*/MBHa. 
Clearly there are two limits: 

• If Mbh ~ TO* there is no loss cone. 

• If Mbh ^ oo then loss cone is fixed. 

For non-isothermal distributions, this is a hard prob- 
lem. Even for flat density, non-isothermal distributions 
the problem becomes non-trivial [31]. 

Clearly for a real cusp composed of a distinct pop- 
ulation of roughly stellar mass objects, the black hole 
moves. We want to estimate the mean free path, time 
scale for wandering and whether the black holes "carry 
the ciisp with them" allowing an adiabatic response of 
the stars to the black hole displacement. We know that 
non-isothermal cusps in general have vbh = V'Hso-bh, 
rj^l {d Chatterjee et al 2002 [31]). 

We can assume ballistic loss cone exit, ie consider the 
problem in the limit where the black hole moves on some 
ballistic trajectory relative to the centre of mass of the 
stars, changing directions and velocity amplitude only 
at discrete intervals, with the change occurring rapidly 
compared to the duration of ballistic motion. A com- 
pact stellar remnant is generally captured if it comes 
within (lOrs) of the black hole. If the the black hole has 
been displaced by a comparable amount, then a "new 
loss cone" has been entered. The time scale for this to 
occur is: 



tic » Wrs/vViso-BH OC Mjfj 

For a ~ IO^Mq BH embedded in a tr « ICQ km/sec 

cusp, we find tic ^ l/rylO^ seconds, or about 1 year. The 
time scale is longer if the mean-free path is shorter, since 
in that case the black hole has to random walk out of loss 
cone, rather than being displaced ballistically. 

The motion should be ballistic on time scales compa- 
rable to some fraction of the cusp dynamical time, very 
roughly, ^ 10^ ~^ years. 

So, orbits with time scale less than 1 year are carried 
with the BH as it wanders. Which are orbits inside ^ 
lOOAU for our canonical IO^MqSMBH. 

This is ~ 10*rs, which is precisely where the dominant 
influx of compact objects is coming from. 

Therefore BH wandering is significant, but possibly not 
dominant, with a critical mass close to lO^M© - more 
massive SMBH wander slowly and carry cusp with them, 
less massive SMBH leave loss cone ballistically. 

It is critical to calculate ry for non-isothermal cusps in 
order to estimate a robust coalescence rate. This may 
require numerical simulations. 

We also want the RMS displacement {rsH), as well as 
the mean free path Tbh- The black hole is effectively 
fixed to the centre when tbh *C rs- 

B. Future broken symmetries 

• Triaxiality. 

If we can't move the SMBH, can we move the stars? 
If the nucleus is triaxial, J is not conserved for in- 
dividual stars on orbits - box orbits and boxlets 
walk through centre. Close to the SMBH, orbits 
are perturbed Keplerian orbits, precessing ellipses, 
the transition is inside rn, at what radius the per- 
turbed Keplerian limit is appropriate is important. 

• Chaos: 

Scattering of stars by SMBH can lead to dynamical 
chaos where the only integral of motion conserved 
for stars moving in the mean potential of the galaxy 
is the energy, ignoring long time scale relaxation 
processes. 

Ergodic orbits — > spherical potential (cf Gerhard & 
Binney 1985). 

However, triaxial systems can be constructed with 
a central SMBH in which chaos is suppressed and 
triaxiality persist to small radii, well inside rh (see 
HoUey-Bockelmann et al 2002). 

Regular -regular scattering, allows continued cen- 
trophilic orbits, and a persistent flux of stars to 
very small radii on dynamical rather than diffusive 
time scales. This is a hard problem, and in par- 
ticular we don't know what real galaxies do yet, 
the existence of theoretical models with suppressed 
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chaos is not sufficient to show that real galaxies 
can achieve this limit, and chaos may dominate the 
evolution of the centres of real galaxies with super- 
massive black holes. 

We need to determine flux to low J from high E 
(large r) and the time scales for evolution of poten- 
tial shape (cf Zhao et al 2002 [33]). 

• Mergers: 

A significant fraction of galaxies has undergone 
mergers with other galaxies, including low mass 
satellite galaxies. Merger processes provide for 
strongly non-axisymmetric time dependent central 
potentials which leads to refilling of the loss cone 
on dynamical time scales, for many cusp dynami- 
cal times. Most of the integrated rate may be from 
post-merger galaxies rather than quasi-static pas- 
sively evolving galaxies. After the stellar popula- 
tion involved in a galaxy merger has settled into 
a stationary distribution, the SMBH brought to- 
gether by the collision may evolve towards coales- 
cence through dynamical processes. At late stages 
of SMBH-SMBH coalescence is driven by gravi- 
tational radiation, and therefore the inner stellar 
cusp around each of the SMBH coming together 
will see a strongly time dependent potential. This 
may drive significant numbers of stars into the grav- 
itational radiation loss-cone leading to substantial 
enhancement of the rate of coalescence of low mass 
objects during the process of binary merger. 

AGN activity may correlate with galaxy mergers, 
and coalescence may thus be enhanced during AGN 
phases. This may lead to significant fraction of co- 
alescences occurring when the SMBH has a sub- 
stantial gaseous disk, and the torquing of the in- 
ner accretion stream during late stages of in-spiral, 
producing a low amplitude periodic flickering of 
the AGN, which would provide a valuable electro- 
magnetic signature for the counterpart to the grav- 
itational radiation emission. 

• Spin: 

Wilson & Colbert [34] noted that major mergers 
may both lead to near maximal spin black holes, 
and cause abrupt changes in the spin parameter 
of the black holes. Black hole spin will change 
the shape and extent of the loss-cone, but at a 
level small compared to other uncertainties. How- 
ever, the waveform of gravitational radiation emit- 
ted during inspiral is sensitive to the spin of the 
SMBH [35]. 

Two issues arise: one is whether we can predict 
the fraction of near maximal Kerr black holes from 
observed statistics of accreting systems and models 
of major mergers; and, whether the observations of 
spin parameters from inspiraling compact objects 
can test models of SMBH evolution. Both issues 
are largely open and need to be explored further. 



• Star Formation: 

The largest outstanding uncertainty arc the IMF, 
the low mass BH mass function, and whether star 
formation is ongoing in galactic nuclei, and if so, 
with what duty cycle., 

• Pop III: 

If the Pop III IMF is peaked strongly to high mass 
[36,37], then at high z, intermediate mass black 
holes are strong contributors to the coalescence rate 
and observable LISA event rate. 

• Black Hole Mass Function: 

We need the zero mass cut-off masses, frequency 
and distribution of masses for low mass black holes 
as function of the IMF and mctallicity (see Fryer 
this proceedings and references therein). Knowing 
the binarity and initial spin of massive stars and 
their compact remnant descendants would be nice 
too. 

• Ongoing star formation: 

If there is ongoing "normal" nuclear star formation 

in inner 1-2 pc of normal spirals, with as little as 
1% duty cycle, then compact objects are replen- 
ished and high coalescence rate sustained. That is, 
if there is star formation in galactic nuclei, and if 
the mass function of the stars formed is similar to 
that of stars formed in the held, or possibly biased 
towards high mass stars, then the nuclear star for- 
mation is important for generating a new popula- 
tion of compact stars in the nuclear regions. This is 
a particularly important process if low mass black 
holes continue to form in the metal rich gas found 
in galactic centres at late times. 

There is some observational evidence for ongoing 
central star formation in the Milky Way [38-40]. 
An open question is whether there is in situ for- 
mation, or off centre clusters that then sink to 
the centre [41-43]. We need to understand the 
formation channel and IMF of stars formed from 
metal rich gas deep in galactic potentials, in possi- 
bly intense radiation environments, near the central 
black holes. 



III. CONCLUSIONS 

• Coalescence of low mass compact objects with 
SMBH is certain. 

• Rates are uncertain but consistent and optimisti- 
cally high. 

• The outstanding model uncertainties are resolvable 
in principle. 

The coalescence of low mass compact stellar remnants 
with central supermassive black holes in galaxies is one of 
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the major hoped for LISA sources, and detection of such 
events with reasonable S/N can provide strong tests of 
both astrophysical processes and fundamental physics. 

LISA will hopefully test some of our models, some 
model parameters will probably remain free until LISA 
constrains them. 
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